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A Bose-Einstein condensation (BEC) has been observed in magnetic insulators in the last decade. The bosons 
that condensed are magnons, associated with an ordered magnetic phase induced by a magnetic field. We 
review the experiments in the spin-gap compound NiCl2-4SC(NH2)2, in which the formation of BEC occurs 
by applying a magnetic field at low temperatures. This is a contribution to the celebration for the 50th 
anniversary of the Solid State and Low Temperature Laboratory of the University of Sao Paulo, where this 
compound was first magnetically characterized. 
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I. I. INTRODUCTION 

Macroscopic systems governed by quantum mechan- 
ics of interacting particles attract a great deal of in- 
terest. Cold atoms and quantum magnets, whose to- 
tal spin is an integer, have interesting similarities that 
show the common physics of these two seemingly differ- 
ent realizations^ - — . All atoms with an even number of 
neutrons satisfy Bose statistics, which accounts for about 
75 per cent of the atoms in the periodic table. Based 
on the Bose-Einstein statistics a gas of non-interacting 
massive bosons condenses below a certain temperature 
Tbec, m which the Bose-Einstein condensation (BEC) 
occurs. This is a macroscopic quantum phenomenon 
characterized by spontaneous quantum coherence persist- 
ing over macroscopic length and time scales. In dilute 
atomic gases this phenomenon was realized experimen- 
tally for cold atoms. Several quantum spin systems in 
solids, which show a magnetic-field induced transition, 
are expected to also show condensation above or below 
a certain critical field 2,7,12 . Studies have shown that the 
magnetic system can be mapped non-locally onto a set 
of weakly interacting bosons on a lattice. 

The sorting phase can be described as a BEC of bosonic 
quasi-particles, in which the magnetic field acts to pre- 
serve the number of bosons. Therefore, the tuning pa- 
rameter to induce condensation in spin-ordered systems 
is not the temperature, but the magnetic field. For parti- 
cles as well as magnons, a macroscopic number of bosons 
condense into a single quantum state — the state of lowest 
energy. The quantum coherence of Bose-Einstein conden- 
sation dates back to the prediction of Einstein, based on 
Bose's work, in 1924. 

In the diluted BEC the macroscopic wave function 
is directly connected with the microscopic energy lev- 
els, providing a complete description of these phenom- 
ena in terms of the Gross-Pitaevskii equation. The con- 
cept of a coherent macroscopic matter wave in inter- 
acting many-body systems is independent of a detailed 
microscopic understanding of particles. The intricacies 
of the many-body problem with interactions that lead 
to non-separable Hamiltonians are solved by this equa- 
tion, which introduces effective potentials that are sim- 



pler than the original interactions, which in turn renders 
the physical problem more tractable. 

Experimental evidence of the BEC in confined weakly 
interacting gases was produced by E. A. Cornell, W. Ket- 
terle, and C. E. Wieman in 1995, leading to a Nobel Prize 
in 2001. That BEC would occur in quantum magnets 
was first predicted in 1991, and several reports describ- 
ing real systems were published after 200Q 7 - i 13 i 14 . The 
possibility of experimental investigation of the BEC in 
quantum magnets has led to a deeper understanding of 
the ground states in strongly correlated systems. Re- 
cent studies concerning the BEC of bosons in magnetic 
systems ranged from a 3D weakly coupled 5=1/2 dimer 
compound, TlCuCl3 T (the first system in which an iden- 
tification was attempted with BEC), a quasi-2D spin sys- 
tem, BaCuSi20^, to a quasi-lD 5=1 chain with single- 
ion anisotropy, NiCl 2 -4SC(NH 2 ) 2 (DTN)i±. In addition 
to these three compounds, other gapped spin systems 
share the same physical picture of the BEC, most of 
which comprise Cu dimers. 

The focus of this article is to review systematic con- 
cepts and experimental realizations of a transition to a 
BEC phase in DTN, where bosons are quanta of mag- 
netic excitations in a magnetically ordered ensemble of 
magnetic moments. Although there is no consensus on 
a specific name for the magnetic excitations, we will use 
the word magnon for the excitations in insulating sys- 
tems generated by the application of a magnetic field. 
This denomination follows the nomenclature used by the 
first papers on BEC in magnetic insulators, since elemen- 
tary excitations in antiferromagnetic systems are called 
magnons 2 -^. 

We hope that the experimental evidence for BEC in 
DTN will stimulate new research in this fascinating quan- 
tum phenomenon. 

The outline of this paper is as follows. Section 2 re- 
views the concepts of the ideal BEC at finite temperature 
and the nature of BEC when magnons are induced by the 
application of a magnetic field. Section 3 presents the 
magnetic properties of DTN, introducing the Hartree- 
Fock- Popov analysis of the hard-core boson Hamiltonian. 
Section 4 presents experimental works that lead to the 
interpretation of the induced phase diagram as a BEC. 
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a) Classical behavior 
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b) Wavelike behavior 




FIG. 1. (a) Classical representation of a solid, a is the average 
distance between the particles, (b) In the wavelike represen- 
tation the position of the particles is given by the thermal 
de Broglie wavelength A. Decreasing the temperatures, the 
wave-packet becomes correlated for A ~ a. For particles in 
the same quantum state, such as bosons at very low tem- 
peratures, the system is represented by a single macroscopic 
wave- function. 



quantum statistics. The density distribution of the con- 
densate is represented by a single macroscopic wave func- 
tion with a well-defined amplitude and phase, which are 
the order parameters of the system, and a Bose Einstein 
condensate is formed. Figure 1 shows illustrates this. 

The Bose-Einstein condensation is associated with 
bosonic particles at low temperatures that condensate 
into a quantum state at a temperature Tbec- At low 
temperatures, the density of bosons with mass m may be 
expressed by the power law dependence equation within 
the BEC picture as 



P = C(3/2) 



/ mk B \ 



3/2 



p 3/2 
BEC-> 



(i) 



where ( is the Riemann zeta function, and ((3/2) = 
2.612±^£. 

For diluted bosons with small interaction, meaning 
that the two-particle scattering length x is much shorter 
than the interatomic distance, the observed Tbec agrees 
well with the theoretical expectation for free particles. 
For less diluted systems, the interaction between bosons 
gives rise to a small shift of TbecA^' 



^(T B ec)/Tbec = cp 1/3 x. 



(2) 



Various authors have attempted to determine the value 
of c, which was found to be in the range 1 — 2, although 
all authors agree on the functional form. 



2.1. Candidates for BEC 



Section 5 summarizes the results obtained for DTN in 
the framework of the BEC of magnons. 



II. II. BOSE-EINSTEIN CONDENSATION 

According to the Heisenberg uncertainly principle 
the position of a particle is smeared out over a dis- 
tance given by the thermal de Broglie wavelength A = 
h/^TrmkBT) 1 / 2 , where ks is the Boltzmann con- 
stant, m is the particle mass, and T is the system 
temperature^^. At room temperature A is typically one 
hundred thousand times smaller than the average spac- 
ing between the particles, a. This means that the matter 
waves of the individual particles are uncorrelated, and 
the system can therefore be described by classical Boltz- 
mann statistics. As the system is cooled, eventually the 
distance between two particles becomes the same as the 
de Broglie wavelength, a « A-iii. The packet-functions 
of adjacent particles overlap, causing the atoms to lose 
their identity, and the behavior of the system becomes 
strongly correlated. The system is now governed by 



There are several candidates for BEC. A first attempt 
to experimentally observe BEC was made with superflu- 
idity. When 4 He is cooled to a critical temperature of 
2.17 K, the liquid density drops, and fractions of the liq- 
uid become superfluid. The strong interaction of 4 He in 
liquid form prevents perfect condensation. The forma- 
tion of Cooper pairs of electrons, which act like bosons, 
is also evidence of the condensation effect, which pro- 
duces super conduct ivityi^. However, in these systems, 
a macroscopic wave function provides only a phenomeno- 
logical description of the superfluid state. 

By contrast, the macroscopic wave function of pure 
condensates is directly connected with the microscopic 
degrees of freedom, providing a complete and quantita- 
tive description of both static and dynamic phenomena. 
This description should cover systems showing a true 
BEC phase. Some systems were cited, such as excitons- 
polaritons, which are light-mass quasi-particles usually 
produced by optical excitations in semiconductors, and 
photons confined in a cavit y 1 i 11 i 15 ~ —. However, excitons 
have a very short lifetime and the photons are absorbed 
by the walls. Because the BEC requires a thermalization 
time to condense in which the number of particles is con- 
served, their characteristics almost prevent these systems 
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from becoming BECs. All these systems are believed to 
be defective condensates. 

Recently it has been reported that a gas of magnons 
at room temperature can be continuously overpopulating 
the lowest energy level, which shows characteristics of 
BEC, when it is driven by microwave pumping . 18 ! 19 

On the other hand, diluted magnons in equilibrium can 
be pure condensates at low temperatures, just like cold 
gas. 



B. 2.2. BEC of Magnons 

Magnetic insulators with magnetic ions are character- 
ized by a gap separating the singlet ground state from the 
lowest-energy excited state. If the gap is closed by the 
Zeeman effect, the resulting quasiparticles can undergo 
the Bose-Einstein condensation transition. 

The BEC phase in quantum magnets is an XY antifer- 
romagnetic state in which the spins spontaneously choose 
a particular orientation in the XY plane. The number of 
bosons is proportional to the net magnetization, i. e., 
the S z component of the spins. On the other hand, the 
number of condensed bosons is proportional to the an- 
tiferromagnetically ordered S x component of the spins. 
Thus the XY antiferromagnetic order parameter maps 
directly onto the order parameter for Bose-Einstein con- 
densation. It has a magnitude (size of the spins in the 
spin language, or number of bosons in the boson lan- 
guage) as well as a phase (direction of the spins). In 
the particle language, the boson number can be tuned 
by the magnetic field. For this to occur, the system must 
be an XY magnet, i. e., must have uniaxial symmetry 
in the plane perpendicular to the applied magnetic field. 
Therefore only certain quantum magnets can be treated 
as BECs. We note that in real systems there will al- 
ways be terms that break the uniaxial symmetry, such 
as dipole-dipole interactions and spin-orbit interactions. 
Consequently the theory is valid only when these terms 
are significantly smaller than the temperature. At very 
low temperatures any real quantum magnet will cross 
over to the Ising universality class. 

Bosons can be achieved by creating and maintain- 
ing a large number of quasi-particles in a system. One 
special case is the creation of magnons at the induced 
magnetic-ordering phase. Figure 2 shows the phase di- 
agram in the temperature-versus field plane of an in- 
duced antiferromagnetic ordering in a field whose mag- 
netic ions are strongly correlated. Magnons can be cre- 
ated and, under special conditions, condensed into the 
ordered phase 2 -^^. 

A procedure for studying the properties of the phase di- 
agram is to understand the transitions that appear when 
the boundaries separating the phases are crossedZ»22r— . 
When a border is crossed vertically to enter the ordered 
phase a thermal phase transition results, driven by ther- 
mal fluctuations. The magnetic and thermal energies 
compete. At temperatures near zero, by contrast, there 
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FIG. 2. The phase diagram from the normal to the BEC state 
can be studied via thermodynamic properties such as magne- 
tization, temperature and specific heat, which are measured 
when the phase boundary is crossed by sweeping the field or 
the temperature, as indicated by the arrows. 



are no thermal fluctuations and a quantum transition re- 
sults, driven by quantum fluctuations. The magnetic en- 
ergy and the zero-point energies compete. This second- 
order quantum transition to the BEC phase is free of 
interactions. However, including even week interactions 
into the description of the condensation has proven to be 
a great challenge, which was overcome by introducting 
effective interactions ^ 10 ! 20 ! 21 . These interactions, well- 
approximated by a contact pseudo-potential, are cen- 
tral to arriving at the Gross-Pitaevskii equation. Under- 
standing the behavior of the phase transition near zero 
temperature is understanding the quantum behavior of 
the system. The transition to the BEC state is charac- 
terized by specific characteristics of the thermodynamic 
properties. 

We studied BEC in magnons because condensation oc- 
curs at temperatures that are much more accessible than 
the condensation of particles. 

Since according to (1) Tbec is proportional to p 2//3 /m, 
we would like to maximize p to observe the BEC at eas- 
ily accessible temperatures. However, that equation as- 
sumes that the bosons are free particles and is therefore 
only valid in the limit of low density. We can overcome 
this problem by considering that Tbec is also depen- 
dent on the boson mass m. By minimizing m, we can 
observe the BEC at relatively high temperatures with- 
out sacrificing the condition of low boson density. In 
quantum magnets, the boson mass is sufficiently small, 
so that the BEC can be realized at temperatures of a few 
Kelvin, in contrast to cold atoms, in which temperatures 
in the nano-Kelvin range are needed to observe the BEC 
in magnetic optical trap o 19 i 21 . 



III. III. MAGNETIC PROPERTIES OF DTN-BEC 

In addition to the properties of NiCl2-4SC(NH2)2, one 
needs well-prepared samples for studies under laboratory 
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FIG. 3. a) Sample of DTN grown in a water solution, b) 
Schematic crystal structure of the compound. The positions 
of the Ni and CI atoms are shown. 



conditions. Single crystals of dichloro-tetra/czsthiourea- 
nickel (DTN) were grown by dissolving a large excess 
of nickel chloride in a saturated solution of thiourea^. 
The crystals grow after about two weeks in a solution 
maintained at a temperature of 35 degrees centigrade. 
Yellow-brown single crystals are tetragonal prisms along 
the c-axis (Fig. 3a). The space group of DTN is 74, 
with two molecules in the unit cell, whose dimensions 
are a = b = 9.558 A and c = 8.981 A. The S = 1 Ni 
atoms form two interpenetrating tetragonal lattices. The 
thioureas are bounded to the nickel by means of the sulfur 
atoms in a square array. The spin must be in an environ- 
ment of approximate uniaxial symmetry (Fig. 3b). The 
metal is in a noncentrosymmetric environment of four- 
fold symmetry. All molecules are oriented in the same 
way and held together by Cl-N hydrogen bonds. The 
Cl-Ni-Cl axes all lie parallel to the crystal c-axis. The 
samples were always chemically analyzed for their main 
constituents. Polarized light allowed us to verify the ex- 
cellent quality of the samples and to check the orienta- 
tion of the crystallographic axes for the measurements, 
as determined by X-ray. The sample masses in our ex- 
periments ranged from a few milligrams to 2 g. 

Initially, to provide magnetic characterization, the 
magnetic susceptibitilies of a single crystal were mea- 
sured in the directions parallel and perpendicular to the 
c-axis, and temperatures ranging from 0.4 to 300 K. 

From these data the parameters applicable to the fol- 
lowing 5=1 spin-Hamiltonian with single-ion anisotropy 
were determined: 

H = J - S J • S J+^ + J2i D ( S D 2 - WbHS?}. (3) 

The obtained parameters are g= 2.30 and D « 9 K. 
Molecular field correction gives the exchange interaction 
parameter J « — 2.2 K (antiferromagnetic). The high 
D value and the anisotropy in the exchange suggest a 



quasi ID behavior. These data were used to anticipate an 
induced ordered phase at low temperature in an applied 
magnetic field below « 11 T. 



A. 3.1. Magnetic Induced Phase Diagram - BEC 

The field-induced phase by level crossing is permitted 
in compounds that follow stringent requirements: in ad- 
dition to having a multiplet above the non-magnetic sin- 
glet, the metal ion must reside in a site of uniaxial sym- 
metry, so that the magnetic field can be applied simulta- 
neously parallel to the anisotropic z axis for all magnetic 
ions. The parameters for DTN and its crystal symme- 
try then yield magnetic properties that can be studied 
within this picture. 

In DTN the high single-ion anisotropy (energy gap) is 
responsible for the splitting of the S = 1 spin triplet of 
the Ni 2+ in a S z = ground state and an S z = ±1 ex- 
cited doublet considerably higher in energy. The excited 
doublet can hop to a neighboring doublet via the trans- 
verse component of the magnetic exchange interaction J, 
resulting in delocalized states analogous to the delocal- 
ized electronic states in crystals. At zero field the small 
exchange interaction between the spins is not enough to 
induce long-range ordering, even at zero temperature. 
The application of a magnetic field H leads to a splitting 
of the doublet, which brings down the component par- 
allel to the field linearly with H until the gap between 
the excited and ground states closes at H = H c . At this 
point magnetic excitations, i. e., bosons, are created (Fig. 
4a). 

The more realistic view of the energy levels that takes 
the interaction into account is shown in Fig. 4b. The 
exchange interaction produces a dispersion of the energy 
levels, and there is now a field range, from H c \ to i^ C 2, 
where the S z = +1 level becomes degenerate with the 
ground state and bosons are created. This gapless phase 
extends throughout this field range. The coupling to the 
magnetic field, which is proportional to the magnetiza- 
tion, controls the boson density. The BEC that appears 
in this intermediate field range corresponds to the coher- 
ent superposition of the S z = +1 doublet and the S = 
singlet at each sit o 14 i 22 ~ — . 

The transition to the magnetic phase can be ade- 
quately described as a condensation of bosons. Three are 
the basic conditions for BEC induced by a magnetic field: 
1) the excited state becoming degenerate with the ground 
state, 2) strong correlation inducing 3D long-range or- 
dering, and 3) a uniaxial symmetry leading to number 
conservation of bosons. 



B. 3.2. Hard-core Boson Model 

The first attempt to interpret a field-induced phase or- 
dering in a BEC was made by Nikuni 7 . In his work the 
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FIG. 4. Schematic diagram (T=0) of the levels of energy, a) 
The lowest energy levels of Ni ++ in an axial crystalline field as 
a function of the external magnetic field, b) The broad band 
indicates dispersion of the lowest excited level due to exchange 
coupling. Condensation occurs in the range between the fields 
Hd and H C 2, in which the magnetization M increases almost 
linearly up to the saturation. 



low- field and low-temperatures regions, where this level 
does not contribute. An analysis taking this level into 
account is discussed in Section 4.1 7 14 i 22 . 

The physics of the magnons depends on the relative 
strengths of the repulsive interaction and kinetic en- 
ergy. When the interaction is attractive, vo < 0, the 
system collapses. If the repulsion interaction dominates, 
the system evolves into a state where the bosons form 
a large lattice that gives rise to a magnetization plateau 
in the M vs. H curve. If the kinetic term dominates, 
the system, which corresponds to a Mott insulator in the 
low-field region, undergoes a condensation at the first 
critical field H& . Above this field the density of bosons 
increases proportionally to the longitudinal component of 
the magnetization M up to the saturation value, M sat: 
at the second critical field H C 2. The result is an an- 
tiferromagnetic state characterized by the presence of a 
staggered transverse moment, M stai perpendicular to the 
applied magnetic field, forming a BEC state. The abso- 
lute value and relative angle of M sta are related to the 
amplitude and the phase of the condensate wave func- 
tion, respectively. In this phase the condensate shows a 
macroscopic order parameter that spontaneously breaks 
global phase symmetry by the emerging antiferromag- 
netic longe-range order. In these conditions, there must 
exist a gapless mode called Nambu-Goldstone in the new 
phase, which is necessary to keep the original invar iance 
of the Hamiltonian. At H C 2 all sites are occupied by 
bosons and the system enters a second Mott insulating 
phased. 



induced phase diagram of TICUCI3 was studied by ap- 
plying the Hartree-Fock-Popov (HFP) mean field anal- 
ysis to the hard-core boson Hamiltonian. For a diluted 
gas of particles the interatomic interaction is sufficiently 
weak, therefore the mean field Gross-Pitaevskii theory is 
a logical tool to study this system. In this model the 
particle scattering is dominated by the two-body contact 
interactions, which are described by the s- wave scatter- 
ing length. The physical implication of this condition is 
that it is highly improbable for more particles to interact 
with each other simultaneously. In hard-sphere gas, the 
scattering length is equal to the diameter of the atoms. In 
the low-density limit only the two particles interactions 
are important. The low-energy Hamiltonian (3) can be 
transformed into boson language via identifying £ + = a) , 
where a) is the boson creation operator, in 

W = ( 6 - V))4 a k + y Yl a Uq a k-q a k a k (4) 

k 

where e = hk 2 /2m is the kinetic energy, vo is the short- 
range repulsion between two magnons that occupy the 
same site and V is the unit-cell volume of the sample. 

In this hard-core model applied to the Ni 2+ ions the 
density of the magnetization along the field direction, S z 
= +1, is mapped onto the boson density. Here the S z = 
— 1 component is neglected, because we are interested in 



C. 3.3. Boson Number Conservation 

At the induced phase the Ni magnetic spins do not 
have any preferred orientation in the plane perpendicu- 
lar to the magnetic field, and the antiferromagnetic order 
is XY-like. The tetragonal crystal structure provides a 
uniaxial crystal- field symmetry about the direction of the 
applied field. It is this symmetry that enforces number 
conservation among the bosons. This forces the bosons 
to remain in the system and macroscopically occupy the 
ground state at low temperatures. Boson- number con- 
servation is a key condition that separates bosonic sys- 
tems that condense from those that do not. Formally, the 
phenomenon of BEC requires the conservation of particle 
number. The boson number must be set by some exter- 
nal constraint or else bosons will be excitations of the 
system and vanish as the temperature is lowered to zero, 
as is the case, e.g., for phonons. In the Hamiltonian (4) 
every creation operator a) is multiplied by a destruction 
operator a. If the Hamiltonian is rotated by an angle (f) 
in the plane perpendicular to the field, then a) — >• a^e 1 ^ 
and a -> ae~^ : such that (a t e^)(ae"^) = a+a. Be- 
cause this Hamiltonian is independent of the number 
of bosons is conserved. Thus, the uniaxial symmetry of 
the Hamiltonian creates a number conservation law for 
the bosons 2 ^. In real situations, the square lattice of the 
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FIG. 5. Effective chemical potential as a function of tempera- 
ture for fermions, bosons and classical particles. For decreas- 
ing temperatures, \i for bosons reaches the zero limit at the 
temperature T cr iticai- As T becomes lower than T cr iticai, M is 
pinned at zero^i 



crystal can introduce a small anisotropy in the ab plane 
along with dipole-dipole interactions or Dzyaloshinskii- 
Moriya (DM) interactions. These effects generally occur 
at lower energy scales and can therefore be neglected at 
the temperatures at which DTN was studied. 21 However, 
experimental studies should be conducted to discard ef- 
fects of this possible anisotropy in the formation of BECs. 

D. 3.4. Low-temperature Limit and Critical Behavior 

The self-consistent HFP approximation leads to the 
following effective Hamiltonian for the magnon system: 

U = ^2(e- [i + v n) a+a fc , (5) 

k 

where n is the density of bosons per magnetic ion, given 
by n = M/M sat . 

The chemical potential 

fi = g^ B (H-H c ) (6) 

determining the number of bosons in the ground state 
characterizes an additional direct contribution to the 
magnon energy from the external field H. 

Bosons condense when the renormalized effective 
chemical potential is zero^i^ (Fig. 5): 

/jLeff = fi - 2v n = 0; (7) 

the factor of two comes from exchange. 

At this "point of condensation", where the gap is 
closed, the applied magnetic field is responsible for driv- 
ing the system to the BEC phase in which bosons are cre- 
ated at (iif c i, T c i). At H c i the magnetization acquires the 



value M C (T) and n c (T) = M C (T)/M sat gives the number 
of bosons created at the transition. Using (6) and (7), 
the low- field boundary of the phase diagram H vs. T can 
be written 

H C (T) - H c (0) = [2v / gf i B ]n c (T). (8) 

On the other hand, the density of created bosons is 
given by Bose-Einstein statistics 

p{T) = {l/V)Y J [V[ex P {e k /T)-l)] (9) 

k 

That is, at the low-temperature limit^i^ 2 - and letting 
the chemical potential be /i e //: 

HT)a<m (^f\ (10) 

since the density of bosons per magnetic ion is n c (T)= 
p(T)V where V is the unit cell volume, (8) becomes 

H C (T)-H C (Q) « [2V5Mb]C(3/2) ( ! ^ff V. (11) 

E. 3.5. Signature of the BEC 

Based on this analysis, the phase boundary can be in- 
terpreted as a BEC of magnons. At low temperatures 
only the lowest excited energy level is considered, so that 
the interaction can be replaced by the constant vq. In 
this case, introducing a hard-core constraint assuming 
an infinite on-site repulsion, only one boson occupies 
each site. Certain properties of the induced magnetic 
phase in TICUCI3 can be explained by the concept of a 
BEC. Plotted as a function of the field, the magnetiza- 
tion shows the cusplike dip at the transition predicted 
by HFP theory ^ 20 1 26 . However, disagreements are found 
with other experiments, wich were later attributed to the 
lack of uniaxial symmetry. 

One of the few compounds entering the BEC phase 
under magnetic fields sufficiently low to be reached using 
standard superconducting magnets is DTN. Most other 
candidates require pulsed fields and yield much less ac- 
curate experimental results. To condensate at T ^ 0, 
magnons should be able to propagate in three dimen- 
sions. For BEC to occur in quasi- ID system, such as 
DTN, there must be weak interaction between the ID 
chains, which permits full mobility of magnons in 3D. 

One of the experimental features of the BEC of 
magnons can be found by analyzing the experimental 
data at the boundary H c i(T), where part of the exci- 
tations condense. 

The signature of the BEC universality class at low tem- 
perature is given by the power-la w 3 i 7 i 12 i 20 i 21 i 25 i 26 

H cl (T) - H cl (0) * AT*, (12) 
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where A is fully determined by the coefficient of (11) as 
a function of vq and the mass of magnons. 

DTN behaves in 3D below 1.2 K as a BEC, which is 
where the long-range order sets in, although the universal 
behavior with <fi — 3/2 is observed only below T < 0.2 K. 
This is because the = 3/2 power-law behavior is a low- 
temperature approximation to the full boson distribution 
function. 
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IV. IV. EXPERIMENTAL DETERMINATION OF THE 
BEC PHASE 

A wealth of useful information about the nature of the 
transition to the ordered phase can be obtained from 
several thermodynamics properties in the temperature 
vs. field plane (Fig. 2). The first thermodynamic prop- 
erty used to study the transitions is the magnetization, 
one of the most fundamental quantities observed in ex- 
periments, followed by specific heat and magnetocaloric 
effect. More properties such as magnetostriction, ther- 
mal transport, EPR, inelastic neutron scattering, and 
other techniques are used to identify the induced mag- 
netic phase as BEC. 



A. 4.1. Magnetization 

The boundary of the induced phase diagram of DTN 
can be experimentally determined by curves of magne- 
tization vs. temperature at fixed iJ, with a cusp-like 
minimum at T c , or by H c at the inflexion in scans of 
magnetization at fixed temperature T. These critical 
points originate from the competition of two behaviors 
of the magnetization, which correspond to the number of 
magnons: above T c , magnetization is determined by the 
number of thermally excited magnons, which decreases 
with decreasing T; below T c the number of condensed 
magnons increases with decreasing temperature, over- 
coming the decrease of thermally excited magnons. As a 
result of this competition on the origin of the magnons, 
M has a minimum at T c , with the value M c (Fig. 6a). 
For the same reason, when the phase boundary is crossed 
by sweeping H at a fixed temperature, the condensation 
occurs at the inflexion field H c (Fig. 6b). It has been 
proved that at H C (T C ) a discontinuity of the magnetiza- 
tion is provided by the sudden vanishing of the chemical 
potential /^^^. 

The magnetization was measured using a vibrating 
sample magnetometer, VSM, in temperatures above 0.4 
K at the High Magnetic Field Laboratory facility of the 
University of Sao Paulo. The sample was immersed in a 
3 He bath, and the temperature was measured by vapor 
pressure and carbon-glass or cernox thermometer. The 
magnetic field was supplied by a NbaSn superconductor. 
Measurement at T = 0.019 K was made using the gradi- 
ent field method in the dilution refrigerator. The transi- 
tions were determined in the magnetization versus field 



FIG. 6. a) Curves of magnetization taken at several fields. 
The minimum indicates the transition temperature T c . b) 
Representative curve of dM/dH. The sharp bending indicates 
the transition field in- 



curves by locating anomalies in the second derivative of 
the magnetization 33 . 



B. 4.2. Specific Heat 

All measurements were made on crystals with the ex- 
ternal field H provided by the 17-20 T superconducting 
magnet in a dilution refrigerator system at the National 
High Magnetic Field Laboratory in Los Alamos. The 
magnetic field was always applied along the tetragonal 
c-axis of the sample. Specific heat was measured with 
two methods, as shown in Fig. 7. 

Measurements at quasi-constant temperature were 
made using a developed ac technique, Fig. 7a. A 
small sample, ~ 1 mg, was mounted onto a home-made 
calorimeter consisting of a Si platform furbished with 
an amorphous-metal thin film heater and a bare Cer- 
nox thermometer. The calorimeter was attached to the 
temperature regulation block and was held in a cryogenic 
vacuum system inside a refrigerator. The sample heater 
is driven by one lock-in amplifier running at frequency 
/, and the temperature is detected with a second lock-in 
locked to the frequency of 2f. In both cases a peak is ob- 
served at the transitions. The observed anomalies in C p 
map the known phase boundaries in DTN quite well. In 
the C p versus H curves a double peak is apparent, with 
the peak near H C 2 far larger than the peak near H&. A 
comprehensive theory of the specific heat has been de- 
veloped. Quantum fluctuations present near iJ c i, but 
not near i^ C 2, reduce the kinetic energy of the magnons, 
causing a strong mass renormalization near the two tran- 
sitions, m{H2) / m{Hi) ~ 3. This effect, closely described 
by analytical and quantum Monte Carlo calculations, ex- 
plains similar asymmetries observed in other properties 
of DTN, such as magnetization, electron spin resonance, 
magnetostriction, and thermal conductivity 3 ^. 

At constant field the method used was the quasiadi- 
abatic heat pulse relaxation technique (Fig 7b). The 
sample was mounted onto a sapphire plate and its 
temperature was monitored by a Ru02 field-calibrated 
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FIG. 7. a) Specific heat as a function of the magnetic field 
for fixed temperatures of 0.40 K ( triangles) and 0.75 K ( 
squares), b) Specific heat in an applied field of 10 T. 
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FIG. 8. a) Magnetocaloric effect determined by monitoring 
T while sweeping H up and down, b) derivative dT/dH for 
several temperatures, where the transition is identified as the 
peak. 



thermometer 14 . 



C. 4.3. Magnetocaloric Effect 

This effect was measured by sweeping the field up and 
down while monitoring the sample temperature with the 
bath temperature held fixed. The device was the same 
that was described in measurements of C p by the relax- 
ation technique 14 . In the magnetocaloric effect, whose 
data we show in Fig. 8a, heating is observed as the mag- 
netic field is swept through the boundary of the phase 
diagram, surrounded by regions of cooling before and af- 
ter the transition. The peak in the first derivative of 
T(H), Fig. 8b, corresponding to the maximum heating 
of the sample, was identified as the phase transition. 
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o Experiment : C V (T,H), AT(H), M(H,T) 
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H c (T)-H c (0)*>AT* 



Scaling 



FIG. 9. Magnetic phase diagram determined from specific 
heat, magnetocaloric effect, and quantum Monte Carlo sim- 
ulations (QMC). The magnetization vs. field measured at 19 
mK and calculated from QMC simulations are overlaid on the 
phase diagram 14 . 



different. In the field-induced quantum critical point, the 
magnetic order is suppressed by both the phase and the 
amplitude fluctuations. Figure 9 shows the constructed 
phase diagram. 

The measured thermodynamic properties of DTN 
C P (H,T), M(H,T), AT(if), used to characterize the 
transitions, were analyzed in regions near the phase 
boundaries, the critical regions, to confirm that this is 
a phase transition to a BEC. Very few quantum magnet 
candidates to BEC have been studied at temperatures 
in which the power-law dependence can be confirmed. 
If the temperature range in which the power-law fit is 
performed is far from zero temperature, it is difficult to 
accurately identify the power-law behavior with correct 
parameters. For example, the first determination of the 
exponent in DTN yielded <f> ~ 2, very different from the 
expected valued. This problem can be circumvented by 
the windowing method, in which the intercept and the ex- 
ponent are determined independently by performing fits 
over different temperature ranges and extrapolating the 
values to zero temperature. The results show that the 
field-temperature phase boundary approaches a power 
law near the quantum critical point, with an exponent 
value (j) « 1.5, as expected for 3D BEG^. 



D. 4.4. Phase Diagram 

The induced ordered phase in the field range H& < 
H < H C 2 extends to a finite temperature T C (H) in which 
the maximum temperature for DTN is T max = 1.2 K. The 
thermal and quantum phase transitions are qualitatively 



Using the theoretical prediction obtained for the scal- 
ing form of the field- and temperature-dependent magne- 
tization close to H C 2(T), an analysis of the magnetization 
shows a good agreement between the theoretical and ex- 
perimental results. This clearly shows that the transition 
at H c2 (T) is a BEC of magnons^. 
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E. 4.5. Neutron Scattering 

A microscopic understanding of the ground state was 
reached through inelastic neutron measurements. Low- 
energy magnetic excitations were performed using the 
cold neutron triple-axis spectrometer at the Paul Scher- 
rer Institute, Switzerland. Using an Oxford instruments 
dilution insert in a VARIOX cryostat, a sample of three 
co-aligned deuterated single crystals of DTN with a com- 
bined mass of 3 g were cooled to T = 80 mK. The analy- 
sis of the resulting dispersion of the magnetic excitation 
was adjusted to the Hamiltonian with interaction. The 
observed scattering was analyzed using a single-mode 
cross section with dispersion obtained from a general- 
ized spin- wave approach for the disordered phase. The 
obtained parameters agree with a model of spins of Ni 
ions strongly coupled along the tetragonal axis, but only 
weakly perpendicular to it, making DTN a weakly cou- 
pled chain system with an anisotropy higher than the 
exchange interactions^. 



F. 4.6. Electron Spin Resonance 

We availed ourselves of the high-field approach exact 
theoretical expressions for the spin-polarized phase. In- 
vestigating magnon excitations in fields up to 25 T al- 
lowed us to obtain a reliable set of spin-Hamiltonian 
parameters. Part of this work was performed at the 
National High Magnetic Field Laboratory, Tallahassee. 
The parameters agree very well with those obtained from 
fitting the experimentally induced phase boundary and 
low-temperature magnetization of DTN with results of 
quantum Monte Carlo simulation. The parameters were 
used to calculate the frequency-field dependence of two 
magnon bound-state excitations predicted by theory, and 
observed in DTN for the first time2Sr^£. 



G. 4.7. Four Sublattice Model 

We present the low-energy excitation spectrum in the 
magnetically ordered phase obtained from spin resonance 
measurements down to 0.45 K. The EPR measurements 
were made at the Kapitza Institute using a transmission- 
type spectrometer equipped with a cylindrical multimode 
resonator and a 3 He cryostat. The observed modes can 
be interpreted within a four-sublattice antiferromagnetic 
model with a finite interaction between two tetragonal 
subsystems, VJ ntra — sub 



H. 4.8. Crossover from the ID Fermionic to the 3D 
Bosonic Character in DTN 

Magnetoacoustic studies in the vicinity of the quan- 
tum critical points in DTN show that the behavior of 



the observed properties outside the ordered phase bound- 
ary can be well described by an effective ID fermionic 
model of low- lying spin excitons. This, together with 
previous results showing the bosonic behavior inside the 
boundary, suggests a crossover from the bosonic to the 
fermionic character of the magnetic excitations close to 
the quantum critical points. One fascinating aspect of 
this crossover is continuous tuning from the bosonic to 
the fermionic statistics. This research was conducted in 
Dresden 39 . 



I. 4.9. Critical Behavior at the Low Field Transition 

Interpretating the low-field phase boundary as a tran- 
sition to the BEC phase allows one to obtain information 
about the boson interaction in DTN. 

The critical field H C1 plotted as a function of n c in 
Fig. 10a, follows the linear relation in (8). The inter- 
action strength constant in that equation can be ob- 
tained from the derivative of the fitted curve dH c /dn c = 
9.52 = 2vo/ g\iB = B as vo — 0.61 meV. The zero tem- 
perature transition field is obtained as H c (0) = 2.15 T, in 
close agreement with Ref. [14]. This confirms the behav- 
ior is linear over the entire temperature range. Using the 
value of A w 0.70(T/ if 3 / 2 ), as obtained from the phase 
diagram at low temperature, and B = 9.52, (8) and (12) 
simultaneously yield the derivative dn c /dT 3 / 2 = A/B — 
0.073, which is plotted in Fig. 10b as a straight line. 
The continuity of this line at low temperature, with the 
experimental points at high temperature, is a robust in- 
dication that the theory is sufficiently coherent to allow 
us to combine results from different sets of experimental 
dataM. 



J. 4.10. Boson Number Conservation 

The derivative dM/dH in Fig. 6b shows a sharp kink 
at H c . In a recent report, a Dzyaloshinskii-Moriya (DM) 
interaction in DTN was proposed to explain the extra 
lines in ESR experiments. This anisotropic interaction 
between corner-center coupling spins at sites z, j in the 
body-centered tetragonal lattice of DTN is given by ~ 
d- (Si x Sj), where d is a specific vector coefficient. For 
any possible DM interaction, the vector d would have to 
point along the tetragonal axis. Otherwise a nonuniform 
field distribution would exist in the sample, which would 
induce a staggered moment perpendicular to the applied 
field, which in turn would broaden the magnetization at 
the transition. This configuration characterizes a bro- 
ken axial symmetry, as is observed in TICUCI3. Traces 
of M vs. H are the best for analyzing the configura- 
tion of the vector d in DTN. From the sharp peak in 
d 2 M/dH 2 at the transition, we conclude that any pos- 
sible U(l) symmetry-breaking terms are small at these 
temperatures, and the number of bosons is conserved, as 
required for BEG^l. 
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FIG. 10. a) Critical transition field as a function of critical 
magnetization density n c . The filled squares represent the 
measurments. The dashed line denotes the linear fitting of 
the data. The Quantum Monte Carlo data are shown as red 
circles, b) n c as a function of T 3 ^ 2 . The straight line is the 
plot of dn c /dT 3/2 = A/B = 0.073, obtained from (8) and (12) 
from different experiments. The continuity of this derivative 
with the experimental points is very high 40 



L. 4.12. 
Class 



Direct Measurement of the BEC Universality 



Although the windowing method yields the critical ex- 
ponent at zero temperature by an efficient extrapolation 
procedure, the physics that can appear below the stud- 
ied temperatures can be hidden by this method. This 
was clearly demonstrated for BaCuSi206, in which a re- 
duction in dimensionality of the spin system changed the 
exponent, as was observed at the lowest temperature. 
In DTN new measurements were performed at tempera- 
tures down to 1 mK, which is two orders of magnitude 
below the lowest temperature scale for magnetic coupling 
in this system, = 0.18 K, and the lowest ever used 
to investigate BEC in a quantum magnet. Thus, with 
this new experiment, no extrapolation was needed to de- 
termine the power-law exponent. Measurements of the 
transition were carried out using ac susceptibility. The 
experiments were performed using a PrNis nuclear refrig- 
erator and a 15 T magnet at the National High Magnetic 
Field Laboratory Facility at the University of Florida. 
The sample was immersed in liquid 3 He in a polycarbon- 
ate cell, rather than being glued to a cold finger, and 
thermal contact to the refrigerator was assured via sin- 
tered silver. The temperature was calibrated using a 3 He 
melting-pressure curve thermometer. The device used for 
operation at ultra-low temperatures is the mutual induc- 
tance bridge assembly. The transitions were determined 
by the derivative of the susceptibility curve. The power- 
law temperature dependence of the phase transition line, 
with (j> = 3/2, was found to be consistent with the 3D 
BEC universality class, which confirms the conclusion 
obtained by the windowing method 42 . 



K. 4.11. Magnetostriction 



M. 4.13. Disorder in the Doped System DTN/Br 



Magnetostriction taken at dilution refrigerator temper- 
atures, performed in a titanium dilatometer, shows sig- 
nificant magnetoelastic coupling and a magneto-order- 
induced modification of the lattice parameters in DTN. 
Length changes were monitored capacitatively with a Be- 
Cu spring- mounted titanium tip in a plastic rotator. The 
top-load dilution refrigerator was mounted in a 20 T su- 
perconductor magnet system at the National High Mag- 
netic Field Laboratory in Tallahasse. A simple theory 
relates the magnetostriction to the spin-spin correlation 
function and describes the data remarkably well. From 
this, was possible extract the spatial dependence of the 
magnetic exchange coupling. The measured data agree 
excellently with those obtained from quantum Monte- 
Carlo simulations and with the phase diagram deter- 
mined from C p and magnetostriction 4 ^. 



When disorder is introduced in a magnet with interact- 
ing bosons, the condensate disrupts and interferes with 
phase coherence. The result is the creation of a peculiar 
state, the Bose-glass phase (BG), with only short-range 
correlation. In the presence of disorder the transition to 
BEC is argued to occur only from the BG, and never 
from the Mott insulator. The transitions to BEC should 
be governed by physical parameters characteristic for BG 
— » BEC for the low-field transition and for BEC — >> BG 
for the high-field transition. Experimental and theoret- 
ical studies were realized in the disorder quantum mag- 
net DTN/Br, Ni(Cli_ ;E Br ;E )2-4SC(NH2)2. The remark- 
able agreement between theory and experiment shows a 
well-controlled realization of a disordered Bose fluid with 
a new seemingly universal exponent governing the scaling 
of the critical temperature from BG to BE C 43 i 44 . 
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N. 4.14. Mass of Magnons 

From the linear behavior of n c at low temperatures in 
Fig. 10b, the effective mass of the bosons in DTN, 7715, 
can be calculated with the coefficient dn c /dT 3 / 2 from (8) 
and (11): 



dn c /dT 3/2 = C(3/2) 



rn b k B 

27TH 2 



3/2 



V = 0.073. 



(13) 



The obtained value is ra&=5.4xl0 -25 g, which is ap- 
proximately 1/3 of the proton mass. The value for rrib 
can be compared with that obtained from neutron dis- 
persion data, m=3.7xl0 -25 g 45 i 46 . 



O. 4.15. Quantum Depletion of Magnons 
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FIG. 11. Magnetization density at T— 0.019 K as a funtion 
of the field. The filled squares represent the measured den- 
sities. The open symbols, representing the condensed bosons, 
are obtained from (14) 



Although the repulsive interaction between bosons are 
essential to avoid the collapse of particles, vq depletes 
the number of condensed bosons by the quantity of non- 
condensed bosons. This was described by Nikuni, who 
obtained the density of bosons in the BEC phase from 
the sum of the condensed and the non-condensed densi- 
ties. To determine the fraction of condensed bosons, we 
analyzed the n vs. H curve at T = 0.019 K in Fig. 11, 
which is reasonable approximation to the T — >• limit. 
The relation between the total density of particles n and 
the density of condensed bosons n con d, as an effect of ^o, 
is given at zero temperature by the equality 



n(0) = n cond (0) + ^ ( ^2 



3/2 



The second term on the right-hand side of this equation 
represents the number of condensed bosons that are scat- 
tered out of the ground state due to interaction ^o- The 
non-condensed fraction at T = 0, called quantum deple- 
tion, £n, is caused by quantum fluctuations around the 
true condensate. Applying (14) to the experimental data 
we determined the dependence of the condensed bosons 
with the magnetic field. The open circles in Fig. 11 
show the calculated density of condensed bosons. The 
depletion of condensed bosons just above the critical field 
H c (0) = 2.09 T increases almost linearly up to « 10 per 
cent at the field of 3.0 T 45 . 

Depletion is an important matter in BEC because the 
dissipation effects of non-condensed bosons leads to the 
creation of solitons and vortices in diluted BECs. In 
DTN, an insulating system, this effect has not been ob- 
served. 



exceptionally well with the experiments. Some measure- 
ments also are correctly fitted by analytic expressions. 
The values for the DTN parameters refined via QMC 
and ESR are D = 8.9 K, J c = 2.2 K and J ab = 0.18 K^. 

With an alternative approach using the diagram tech- 
nique, one finds an expression for the magnetic phase 
boundary. Taking into account the coupling constant 
between the two different tetragonal sublattices, the fol- 
lowing set of parameters fit the experimental data very 
well: D = 7.72 K, J c = 1.86 K, J ab = 0.2 K, and 

Mntra-sub « 0.1 K^. 



(14) Q. 4.17. Semi-hard -core Boson Model 



Both TICUCI3 and BaCuSi206 have a strong antiferro- 
magnetic coupling between the two 5=1/2 spins in each 
Cu dimer, forming an 5=0 singlet ground state sepa- 
rated by a gap J from the 5=1 excited triplet. The 
hard-core model describes the experimental symmetry- 
induced phase diagram of these systems very well. In 
DTN this model describes the behavior of the transitions 
only at low temperatures, T < 0.2 K, and low fields, 
H < H c \. An alternative approach can be used to de- 
scribe the asymmetry of the phase diagram at high fields 
mapping the 5 = 1 system onto a gas of semi-hard-core 
bosons. This asymmetry is attributed to the influence of 
the highest energy level and is not expected for the ef- 
fective model with 5=1/2. In this model the maximum 
boson population is limited to two per lattice site. More 
studies focus on the analysis of the upper critical field 
and a self-consistent theor y 48 ! 49 . 



P. 4.16. Theoretical Methods 



R. 4.18. Two-magnon Bound States 



Almost throughout, the measurements are compared 
with quantum Monte Carlo data. The simulations agree 



Theoretical and experimental systematic high-field 
ESR studies of DTN have been presented with special 



12 



emphasis on single- ion two-magnon bound states. To 
clarify some remaining discrepancies between theory and 
experiment, we analyzed the frequency- field dependence 
of the magnetic excitations in this material. In particu- 
lar, a more comprehensive interpretation of the experi- 
mental signature of single-ion two-magnon bound states 
is shown to be fully consistent with theoretical results. 
Moreover, we have clarified the structure of the ESR 
spectrum in the so-called intermediate phase. The exper- 
imental part was performed at the National High Mag- 
netic Field Laboratory, Tallahassee. 50 



V. V. SUMMARY 

We presented an experimental review of the recently 
observed behavior of DTN associated with BEC phenom- 
ena. These experiments provide opportunities for an in- 
troductory study of the BEC of magnons in insulating 
compounds when bosons are induced by an applied mag- 
netic field. The simple magnetic structure of DTN en- 
ables experiments in field and temperature regions that 
are easily accessible in the laboratory. As shown in the 
theoretical outline presented here, the many-body as- 
pects of BEC are reduced to an effective single-particle 
description, where interactions give rise to an additional 
potential proportional to the local particle density. This 
model allows characterizing the condensation of magnons 
at low temperatures through parameters such as critical 
exponents, analysis of the boson density, boson inter- 
action strength, boson mass, depletion, etc. Although 
this hard-core model describes the low-field behavior of 
DTN very well, the whole phase diagram can be exactly 
explained by quantum Monte Carlo simulations with a 
semi-hard-core model. 

The purpose of this review was to describe experiments 
that go beyond the mean-field theory and explore effects 
that modify the transition to the field-induced magnetic 
phase. 

Finally, we note that the research in this field has not 
been limited to condensed matter, but was extended to 
the applications in nuclear physics, astrophysics, and par- 
ticle physics. 
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